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Abstract 

We show that the inclusion of long distance effect, expected in 
strong interactions, through a nonextensive thermodynamical approach 
is able to explain the experimental distribution of the transverse mo- 
mentum of the hadrons with respect to the jet axis (pt) e^e~ — > 
hadron reaction. The observed deviation from the exponential behav- 
ior, predicted by the Boltzmann-Gibbs thermodynamical treatment, 
is automatically recovered by the nonextensive Tsallis statistics used 
here. We fitted the observed pt spectrum in the range of 14 GeV to 
161 GeV and obtained, besides a good fit, the theoretical important 
fact that the temperature becomes independent of the primary energy. 
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1 Introduction 



The global structure of multiple hadropro duct ion in e^e~annihilation, has 
been well understood with remarkable results concerning the total cross sec- 
tion and angular distribution calculations. Perturbative Quantum Electrody- 
namics and perturbative Quantum Chromodynamics provide a good descrip- 
tion of the initial process involving short distances e'^e" — quark- ant iquark 
(qq) and e^e~ quark- ant iquark-gluon (qqg) interactions, main responsible 
for the global characteristics of this kind of process. The quarks produced 
off-shell initialize a cascade of several quarks and gluons with complex inter- 
actions themselves. This stage is in general described in probabilistic terms, 
based on leading-log Quantum Chromodynamics approximation. The out- 
going colored partons are transformed into color singlet hadrons through the 
soft hadronization process, forming a jet of particles traveling approximately 
in the initial parton direction. However, hadronization require the Quantum 
Chromodynamics in the soft regime, where the coupling constant become 
large and there is no easy way to understand it from first principles. 

For more than a decade several different models making use of Monte 
Carlo techniques were developed to describe such processes [||, 0, |^. These 
methods avoid the microscopic complexity by being able to define probabil- 
ity distributions for each stage of the jet evolution. In these approaches, all 
known aspects involving Quantum Chromodynamics were included, the un- 
known aspects like hadronization being modeled phenomenologically. Each 
model has a set of free parameters which enables it to reproduce many distri- 
butions over a wide range of energy. In particular, the transverse momenta 
of the hadrons with respect to the jet direction, main subject of this paper, 
has been accurately parameterized. Nevertheless, there is no fundamental 
understanding of the matter . 

One alternative approach to understand the transverse momentum distri- 
bution of the charged hadrons produced in e"'"e~annihilation is to look at the 
whole process in a macroscopic way, trying to extract information through 
a thermodynamical treatment, not taking into account the microscopy in- 
teractions between partons, governed by the Quantum Chromodynamics []. 
Fermi was the pioneer in this kind of approach, followed by Hagedorn 

^We used the total transverse momentum instead of pj" and because we are looking 
at the global characteristics of the events 
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who provided the most consistent and well succeeded model. Nowadays, the 
thermodynamical formalism is widely used in heavy ions collisions and 
has also been applied in e^e~, proton-proton (pp) and proton-antiproton 
{pp) interactions to determine the multiplicity of hadrons produced, with 
noteworthy results 0, |^ . 



2 Statistical approach 

According to Hagedorn, a high energy collision produces an unlimited 
and undetermined number of different excited hadron fireballs that reach a 
thermodynamical equilibrium. In this relativistic statistical Ansatz, not only 
the number of particles but also the number of available types of particles 
itself grows with energy. The immediate and important consequence of this 
production mechanism is that the temperature is independent of the primary 
energy p. A similar statement was also made, a few years later, by Field 
and Feynman [^. One should expect that this temperature governs the 
transverse momentum (pt) distribution of the outgoing particles with respect 
to the jet axis. 

In a modern Ansatz [0 , one can assume that the thermodynamical equi- 
librium is reached in each jet produced in e~^e~ — > hadrons interaction. We 
consider in this paper that the thermalization begins just after the hadroniza- 
tion process takes place and finishes when the hadrons are no longer inter- 
acting. Since the range of this interaction is of order of a few Fermi, the 
resonances produced, with a typical lifetime of order of 10~^^ s, should decay 
while the particles are still interacting. Thus, the particles produced from 
these decays should also participate in the thermalization process. As most 
charged particles observed in the detector are produced prompt or through 
the decay of hadronic resonances, we could expect that the charged particles 
observed in the detectors have their distributions dictated by the previous 
'thermodynamical' equilibrium. 

It is then important to determine the thermodynamical properties of the 
event, without any kinematical effect caused by the fast relative motion 
between the hadrons. Using the Boltzmann-Gibbs statistics in this large 
grand-canonical ensemble for a relativistic gas, the transverse momentum 
distribution with temperature Tq, can be written as (see Hagedorn p): 
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where pi means the longitudinal momentum, //^ = p^+m^ and m is the mass 
(m ^ pt) @1- After the integration indicated in eq. (0) and for pf/To ^ 1 we 
obtain an exponential asymptotic distribution for the transverse momentum 

~cpi^exp(-|f) . (2) 
a dpt To 

As the temperature should not depend on the center-of-mass energy {Eqm), 
we can expect that eq. (H) presents the same shape for all e~^e~ energies 
spectrum varying only the parameter c, which depends on the average mul- 
tiplicity of the events and consequently on the Ecm- 

The measured transverse momentum distributions of charged hadrons 
with respect to the jet axis, defined as the sphericity axis |12|, for center- 
of-mass energies in the range of 14 GeV to 161 GeV are displayed in figure 
1. One notices the clear experimental deviation, for high values of pt, from 
the exponential behavior predicted by eq. (||) and indicated by a dotted line. 
Notice that we could also get decent fits for low values of data point of pt 
{pt < 1.5 GeV) using the eq. (0), but in this case the best values for the 
parameter Tq would clearly increase with energy, a feature which disagrees 
with Hagedorn's physically reasonable assumption. 



3 Alternative statistical approach 

The Boltzmann-Gibbs statistics, used to derive eq. (|^), is suitable for an 
enormous number of systems, but it is known that it could lead to unphysical 
results for some systems that present long-range interactions (LR) |. In fact, 
LR interaction is the usual interpretation for the intermittency phenomena 



[|1^, large particle density fluctuations |]14[ |15|, present in hadroproduction. 
Long-range (LR) interactions can occur in many situations, particularly in 
processes where the interacting volume is of order of the dimension of the 

^Note that we are using the expression - long range interactions - in the statistical 
mechanical sense, meaning interaction among particles that are far away one of the other. 
This situation include both cases, soft and hard interactions in QCD. 
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hadrons themselves, resulting in a high hadronic density. In this way the 
gradual deviation of the observed distribution from the exponential function, 
shown in figure 1, could be interpreted as a gradual increase of the multiplicity 
and the hadronic density in the thermalization mechanism. 

Therefore it is useful, to look for another statistics that, unlike the Boltzmann- 
Gibbs, could treat these systems. Recently Tsallis [0 proposed a gener- 
alized statistics whose formalism |T6| , [T^ has been able to describe many 
systems having as common features the presence of LR interactions and/or 
non-Markovian processes [jl8|. This formalism is based on the entropic form 
Sq = {1 — Trp'^)/{q — 1) where p is the density matrix (state) of the system 
and g is a real positive parameter. When g — > 1 the above entropy tends to 
the well-known Boltzmann-Gibbs-Shannon (BGS) entropy S = —Trplogp. 
In this sense the Sq entropy extends the BGS entropy to the nonextensive 
cases corresponding to g 7^ 1. In fact, in this case the additivity property of 
the entropy fails because it appears a new term proportional to g — 1. For 
example, if we look the way the composed entropy of two sub-systems A and 
B, Sa+b, is related to the entropy of each of the two sub-systems, 5*^ and 
Sb, when the sub-systems A and B are statistically independents, one gets 
the expression Sa+b = 5'a + + (1 — q)SASB- It is easy here to see that the 
case g = 1 reobtains the standard result (additivity) from Boltzmann-Gibbs 



entropy |T7| , pO| . The extremization of the Sq entropy with adequate con- 
straints lead us to a probability distribution that exhibits a power-law decay 
(instead of the exponential decay presented by the BGS entropy). Examples 
of application of this statistics are available in several papers |T^, whose 
studied systems are, in general, nonextensive. Thus, assuming a thermal 
equilibrium of a relativistic gas (the fireball) that obeys the Tsallis statistics 
we obtain, observing the derivation made by |T^, the following equation for 



the transverse momentum distribution (which replaces eq. (|I])): 
1 da f°° J I. 1 - 9 



cpt / dpi[l - — ^Vp? + /i^] , (3) 
a dpt Jo Tq ^ 

where g is related to the degree of non-extensivity of the system |T^, |T^, pp] ]. 
To, pi, c and p have the same meaning as in the previous equations. The 
Boltzmann-Gibbs limit - eq. (|l]) - is achieved as g ^ 1. The exact expression 
of the transverse momentum distribution, eq. (|^), is written in the appendix. 
With this expression we fit the experimental data leaving g, Tq and c as 
free parameters and the best fits are presented in the figure 1. We can see 
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Experiment 


Energy (GeV) 




To(GeV) 


TASSO 


14 


1.020 ±0.005 


0.130±0.002 


TASSO 


34 


1.1225 ±0.0006 


0.1152 ±0.0004 


DELPHI 


91 


1.1938 ±0.0005 


0.1094±0.0004 


DELPHI 


161 


1.215 ±0.002 


0.110 ±0.002 



Table 1: 



that the asymptotic behaviour of the new distribution, for larger values of pt, 
represents very-well the deviation from the exponential for all center-of-mass 
energies. 

In figure 2 we present q and Tq values as functions of the e+e~center- 
of-mass energy. Their numerical values are shown in table 1. We can see 
that the q parameter increases smoothly with the energy, saturating around 
1.2. On the other hand, the temperature has a pretty constant behaviour, 
To ~ 0.11 GeV, for the higher energy range. We note that this constant 
value for the temperature is not an ad hoc hypothesis; rather, it was ob- 
tained by the fit, indicating that all the energy given to the e^e~system is 
being used to create new fireballs, as predicted by Hagedorn |5|. We note 
that our Tq value is lower than the value obtained in reference 0, and this 
difference is probably due to the fact that, unlike us, the author considered 
the thermal equilibrium just after the hadronization, i.e., before the decay 
of the hadronic resonances. In our process, the interaction volume continue 
to expand, lowering the temperature, until the particles no longer interact 
among themselves. Also, the statistics used in both papers are different. The 
smooth increase of q with the center of mass energy can be understood as the 
expected increase of the influence of the LR interactions among the hadronic 
particles produced in the event. 

The deviation from an exponential behaviour observed in the transverse 
momentum distribution in the e"'"e~interaction, is also observed for hadrons 
produced in other interactions like pp or heavy ions. It indicates a universal 
characteristic of the dynamics of hadron production, probably due to the 
relevance of the LR regime, related to the high hadronic density. 
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4 Conclusion 



The thermodynamical treatment using the TsaUis statistics discussed here 

presented a remarkable agreement with the experimental data in more than 
4 orders of magnitude in the differential cross section - emphasized with the 
predicted constant behavior of the temperature Tq. Therefore, the results of 
this approach suggest an alternative statistical treatment for the transverse 
momentum distribution within the accepted theory of strong interaction. 
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Appendix 



After some manipulation, the solution of the integral in eq. was obtained 
by the software Mathematica 3.0 and the transverse momentum distribution 
can be written as: 



(4r (I - r - 2) * 

p ( _g o _2_ _2 3 (g-l)pf+T \ I 

2-^ 1 V«-l '(/-I' (/-I 2' 2{g-l)pt y 

p f_g o _2_ _2 1 (g-i)pf+T \ , 

2-^11,9-1 '^'g-l'g-l 2' 2{g-l)pt ^ 

r 9 19 g 1 (g-l)Pt+T \N I 

MFT " 2^1 (^pT - i, ^ - 2' 2(q-l)p, JH 

p f _2 1 _g_ _2 L i (g-i)pt+r \xx 

2-^ Hg-l 'g-l'g-l"^2' 2(q-l)pt ) > > 

2F,Regulamzed (-i, f , | - ^, ^gl^gf ) + 

2 (^) p,2Fi/2e^?«/anzed (i, |, | - ^, ^f(^) - 

3 (l - (^) Pt) * 

,F,Regulanzed (|, |, | - ^, ^^^^)))Sec (^)))/ 
(l6/(^ (l + ^Pt) ^ (^) (g - 1)^) 

where 2-^1(0, 7; means the generalized hyper geometric series, 
2FiRegularized{a, z) = 2Fi{a, z)/T{j) and 
T{x) is the Gamma function. 
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Figure 1: Transverse momentum distribution . The distribution --j^ of the 
transverse momentum pt of charged hadrons with respect to jet axis (de- 
fined in these experimental results as the sphericity axis) is sketched for four 
different experiments, whose center of mass energies vary from 14 and 34 
Gev (TASSO) up to 91 and 161 Gev (DELPHI). The Hagedorn predicted 
exponential behaviour is shown by the dotted line. We can see that the devi- 
ation of the exponential behaviour increases when the energy increases. The 
continuous lines are obtained from our eq. (^ and agree very-well with the 
experimental data. The inset shows the transverse momentum distribution 
for small values of p^. 
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Figure 2: Variation of (a) the entropic index (q) and (b) the temperature 
(To) with the center of mass energy. The entropic index q (stars) increases 
shghtly and the temperature (squares) remains almost constant, as predicted 
by Hagedorn, . The center-of-mass energies vary from 14 to 161 Gev. The 
error bars are smaller than the size of the symbols. 
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